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Deriving Boltzmann Equations from Kadanoff— Baym Equations 

in Curved Space— Time 

A. Hohenegger|j A. KartavtsevL] and M. Lindnei[3 
Max-Planck Institute fiir Kernphysik, Saupfercheckweg 1, 69117 Heidelberg, Germany 

To calculate the baryon asymmetry in the baryogenesis via leptogenesis scenario one usually uses 
Boltzmann equations with transition amplitudes computed in vacuum. However, the hot and dense 
medium and, potentially, the expansion of the universe can affect the collision terms and hence the 
generated asymmetry. In this paper we derive the Boltzmann equation in the curved space-time 
from (first-principle) Kadanoff-Baym equations. As one expects from general considerations, the 
derived equations are covariant generalizations of the corresponding equations in Minkowski space- 
time. We find that, after the necessary approximations have been performed, only the left-hand 
side of the Boltzmann equation depends on the space-time metric. The amplitudes in the collision 
term on the right-hand side are independent of the metric, which justifies earlier calculations where 
this has been assumed implicitly. At tree level, the matrix elements coincide with those computed 
in vacuum. However, the loop contributions involve additional integrals over the the distribution 
function. 

Keywords: Kadanoff— Baym equations, Boltzmann equation, curved space— time, expanding universe 



I. INTRODUCTION 

As has been shown by A. Sakharov [l|], the observed 
baryon asymmetry of the universe can be generated dy- 
namically, provided that the following three conditions 
are fulfilled: violation of baryon (or baryon minus lep- 
ton) number; violation of C and CP; and deviation from 
thermal equilibrium. 

The third Sakharov condition raises the question of 
how to describe a quantum system out of thermal equi- 
librium. The usual choice is the Boltzmann equation 
01 S 1J, 0- However, it is known to have several short- 
comings. In particular classical Boltzmann equations ne- 
glect off-shell effects, introduce irreversibility and fea- 
ture spurious constants of motion. A quantum mechani- 
cal generalization of the Boltzmann equation, free of the 
mentioned problems, has been developed by L. Kadanoff 
and G. Baym [6|. Direct numerical computations demon- 
strate that already for simple systems far from thermal 
equilibrium the Kadanoff-Baym and Boltzmann equa- 
tions do lead to quantitatively, and in some cases even 
qualitatively, different results 0, H, i, [13, [IH, [ll ■ Study- 
ing processes responsible for the generation of the asym- 
metry in the framework of the Kadanoff-Baym formalism 
is therefore of considerable scientific interest. 

The application of the Kadanoff-Baym equations to 
the computation of the lepton and baryon asymmetries in 
the leptogenesis scenario p^] has been studied at different 
levels of approximation by several authors [iJi Il3 ^'^'i 
lead to qualitatively new and interesting results. How- 
ever issues related to the rapid expansion of the universe, 
which drives the required deviation from thermal equilib- 



rium, have not been addressed there. The modification of 
the Kadanoff-Baym formalism in curved space-time has 
been considered in [l^ [l3, [H, [3 ; where it was applied 
to a model with quartic self-interactions and a 0{N) 
model, though the dynamics of quantum field theoretical 
models with CP violation remained uninvestigated. 

Our goal is to develop a consistent description of 
leptogenesis in the Kadanoff-Baym and Boltzmann ap- 
proaches and to test approximations commonly made in 
the computation of the lepton and baryon asymmetries. 
In particular, we want to find out how the dense back- 
ground plasma and the curvature of spacetime affect the 
collision terms of processes contributing to the generation 
and washout of the asymmetry, check the applicability of 
the real intermediate state subtraction procedure in the 
case of resonant leptogenesis [20, [2l[ , and investigate the 
time dependence of the CP-violating parameter in the 
expanding universe [l5|. 

Since this is a rather ambitious goal, we first study a 
simple toy model of leptogenesis containing two real and 
one complex scalar fields, which mimic the heavy right- 
handed Major ana neutrinos and leptons respectively [22l |. 
The peculiarities of the calculation, related to the pres- 
ence of a gravitational field, are determined only by trans- 
formation properties of the quantum fields - scalar fields 
in this case. For this reason, in the present paper, we 
use a model of a single real scalar field with quartic self- 
interactions, minimally coupled to gravity, to illustrate 
the main points. That is, we use the Lagrangian 
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which does also have the advantage, that one can com- 
pare the derived equations with their Minkowski space- 
time counterparts [3] and with the results obtained in 
[H, [O, [H, [ll|. The formahsm presented here will be 
used to analyze the toy model of leptogenesis. 

The starting point of our analysis, which is mani- 
festly covariant in every step, is the generating func- 



tional for the (connected) Green's functions. Performing 
a Legendre transformation we get the effective action, 
which we use to derive the Schwinger-Dyson equations 
in Secini These are equivalent to a system of Kadanoff- 
Baym equations for the spectral function and the statisti- 
cal propagator, which we derive in Sec. lIVI Employing a 
first-order gradient expansion and a Wigner transforma- 
tion we are lead to a system of quantum kinetic equations 
which we study in Sec.|Vl Finally, neglecting the Poisson 
brackets and making use of the quasiparticle approxima- 
tion, we obtain the Boltzmann equation in Sec. lVIl 



reads |2E 



• The Kadanoff-Baym equations and the derived 
Boltzmann equation are covariant generalizations 
of their Minkowski-space counterparts. 

• The space-time metric enters its left-hand side in 
the form of the covariant derivative, whereas the 
collision terms on the right-hand side are indepen- 
dent of the metric. 

• At tree-level the collision terms coincide with those 
calculated in vacuum, whereas the loop corrections 
contain integrals over the distribution function. 

• In the loop contributions one can clearly distinguish 
the initial, final and on-shell intermediate states, 
which is not the case in the canonical formalism. 

We discuss these results in more details and draw the 
conclusions in Sec. lVIII 



II. SCHWINGER-DYSON EQUATIONS 

In the derivation of the Schwinger-Dyson equations we 
employ results from [1^, [23, ^M ■ Our starting point is the 
generating functional for Green's functions with local and 
bi-local external scalar sources J{x) and K{x,y), 

Z[J, K] = JS^ip exp[i{S + J(p+ ^^Kif)], (2) 

where the action S is given by the integral of the La- 
grange density over space. The Minkowski space-time 
volume element df^x is replaced in curved space-time by 
the invariant volume element y/—gd'^x, where y/—g is the 
square root of the determinant of the metric: 

S — fyZ—gd^x C . 

In the Friedmann-Robertson- Walker (FRW) universe we 
have ^/—g = a^irf)^ where a is the scale factor and 77 deno- 
tes conformal time. The invariant volume element enters 
also in the scalar products of the sources and the field 

gd xJ{x)ip{x) , (3a) 

'gd'^xy/^d'^y (p{x)K{x, y)(p{y) . (3b) 

The functional integral measure is modified in curved 
space-time as well. For scalar densities of zero weight it 
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The evolution of the quantum system out of thermal 
equilibrium is performed in the Schwinger-Kcldysh for- 
malism [25I, [2al . In this approach the field and the ex- 
ternal sources are defined on the positive and negative 
branches of a closed real-time contour, see Fig.[l] the 
functions^ on the positive branch being independent^ of 
the functions on the negative branch. This applies also 
to the metric tensor, i.e. 5+ 7^ gT^^, in general. 



t 



FIG. 1: Closed real-time path C. 

In realistic models of leptogenesis the contribution of 
the heavy right-handed neutrinos to the energy density of 
the universe is less than 5% and can safely be neglected. 
In other words, leptogenesis takes place in a space-time 
with a metric, whose time development is (in this approx- 
imation) independent of the decays of the right-handed 
neutrinos and determined by the contributions of the ul- 
trarelativistic standard model species. Correspondingly, 
in our analysis of the toy model of leptogenesis, we will 
also neglect the impact of the scalar fields on the expan- 
sion of the universe^. This implies in particular that the 
metric tensor on the positive and negative branches is 
determined only by the external processes, and one can 
set ff/ti/ = ff/Ti/ = 9fi^- To shorten the notation we will 
also suppress the branch indices of the scalar field and 
the sources. 

The existence of the two branches also affects the defi- 
nition of the S function: 6{x, y) is always zero if its argu- 
ments lie on different branches [28| . In curved space-time 
it is further generalized to fulfill the relation 



~gf{y)5^{x,y) = f[x), 



(4) 



where the integration is performed over the closed con- 
tour. The solution to this equation is given by [23| 



5S(x,y) = {-g^) ^5{x,y){-gyy 



(5) 



^ In particular there are two local (J4. and J_) and four bi-local 

{K^j^, K^ , -R'-+ and K ) sources. Analogously, the field 

value on the two branches is denoted by (/3_|_ and ip— respec- 
tively, whereas the two— point function components are denoted 
by G++, G+-, G-+ and G JH. 

^ With the exception of the point t = tmax ■ 

^ A theoretical analysis of the back— reaction of the fields on the 
gravitational field has been performed in II6II . An analysis, with 
very interesting numerical results, of a model with quartic self- 
interactions in the Friedmann— Robertson-Walker universe has 
been carried out in [l9ll . 



The generalized 5 function is used to define functional 
differentiation in curved space-time [29*| 

5(j){y) £™ e 

From the definition ([5]) it follows immediately that 
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The functional derivatives of the generating functional 
for connected Green's functions 



W[J,K] = -i\nZ[J,K] 
with respect to the external sources read 
dW[J,K] 



(8) 



dJ{x) 
dW[J,K] 
dK{x,y) 



H^) , (9a) 

i[G(y,x) + $(a;)a>(y)], (9b) 



where $ denotes expectation value of the field and G 
is the propagator. The effective action is the Legen- 
dre transform of the generating functional for connected 
Green's functions, 

r[$, G] = W[J, K]-J<^- ^tT[KG] - i$i^$ . (10) 

Its functional derivatives with respect to the expectation 
value and the propagator reproduce the external sources: 



ST[G,<i>] 

S^{x) 
(5r[G,$] 



= -J{x)- J^^d^z K{x,z)^{z), (11a) 



hKiy,x). 



(lib) 



SGix,y) 
Next, we shift the field by its expectation value 

ip -^ ip + ^ . 
The action can then be written as a sum of two terms 

5M^5e/[$]+5[(^,$]. (12) 

Sci denotes the classical action, which depends only on 
$, whereas S'[(y9, $] = Soif] + Sint[f,^] contains terms 
quadratic, cubic and quartic in the shifted field ip. The 
free field action can be written in the form 

•^o = i JJV^.,d'x^^yd'y p (z^-i) ip , (13) 

where ^~^ is the zero-order inverse propagator 

^-i(a;, y) = t{a^ + hP) 5<^{x, y), D, = g^.V^V^ . (14) 

Since the integration measure in the path integral is 
translationally invariant, the effective action can be 
rewritten in the form 

r[$,G'] ^ - i\n ^ S!peicp[i{S + Jp + \LpKip)] 

+ Scim-\tr[KG]. (15) 



Now we tentatively write the effective action in the form 
r[$,G] EE Sci[<^] + flndet [G-^] + §tr [^-^G] 



r2[$,G] 



(16) 



det 



defining the functional T2- The third term on the right- 
hand side is defined by 

tr p-^G] = J/y^^d^xy^d^y ^-i(a;, y)G(y, x) , 

whereas the second term on the right-hand side is defined 
by the path integral 

G~^1 

-^ =/^(^exp((/7G~V) • 

Using pTJ) we can find the functional derivatives of F. 
Differentiation of tr [^^^G] with respect to G is straight- 
forward and gives 

^-t,[^-^G]='^-\y,x). (17) 

oG{x,y) 

To calculate the functional derivative of In det [G^^] we 
take into account that in curved space-time 

JV^d'^zG~'^(u,z)G{z,v) = SS{u,v) . (18) 

After some algebra and use of ([T8|) we obtain a result 
analogous to that in Minkowski space-time 

S 



■lndet[G-i] =-G-\y,x). 



(19) 



SGix,y) 

The functional derivative of (fT6)l with respect to G then 
reads 



<5F[G,$] 



\G-'^,x) 



'Hy,x) + 



ST2[G,^] 
6Gix,y) 



SG{x,y) 

= -lK{y,x) (20) 

Solving pp]) with respect to K and substituting it into 
(I16p we can rewrite the effective action in the form 



F2[G,$] = -Jln I S!ipexp 



S- 



sr2 



tr 



,5F2 



G 



5G 
where again S = Sq 



: In det [G ^] + const. 



(21) 



Sint , but now with 5*0 given by 



'S'o = i nV^xd^^V^yd*y f{iG ^)if. 



(22) 
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This implies that iF2 is the sum of all 2PI vacuum di- 
agrams with vertices as given by Cint and internal lines 
representing the complete connected propagators G [30] ■ 
Physical situations correspond to vanishing sources. 
Introducing the self-energy 

.(5F2[G,$] 



n{x,y) 



2i- 



(23) 



SG{y,x) ' 
we can then rewrite (j20p in the form 

G-\x,y)^^-\x,y)-nix,y). (24) 

Thus the above calculation yields the Schwinger-Dyson 
(SD) equation. Let us note that the derived equation has 
exactly the same form as in Minkowski space-time. 



III. 2PI EFFECTIVE ACTION 

The structure of the Schwinger-Dyson equation is de- 
termined only by the particle content of the model (here 
a single real scalar field) and completely independent of 
the particular form of the interaction Lagrangian. The 
latter determines the form of the 2PI effective action. 
The lowest order contribution is due to the two-loop dia- 
gram in Fig. [21 which only takes into account local effects 




FIG. 2: Two-, three-, and four-loop contributions to the 2PI 
effective action. 

and cannot describe thermalization. Thus one usually 
also considers the three-loop diagram, which describes 
2 <-> 2 scattering. In addition we take into account the 
four-loop contribution. As is demonstrated below, in the 
Boltzmann approximation it describes the one-loop cor- 
rection for 2 <-!■ 2 scattering. The resulting expression for 
the effective action is similar to that given in @, @, [H, [3l| : 



zr2[G] = ^zr(")[G], 



(25) 



iT'~^\g] 



.(3)r 



A2 



j^f^^d'^xG'^{x,x), 






iT^2\G] = -^Jy^^d^Xy^yd^yy^J^Z 

xG\y,x)G^{x,z)G^z,y). 

Note, however, the presence of the \/—g factors which 
ensure invariance of the effective action under coordinate 
transformations. 

Using the definition of the self-energy ((23)) and the 
functional differentiation rule in curved space — time we 
obtain 



n(x,y)=^n(")(x,y). 



(26) 



n(2)(x,2;) 



iS^x,y)-G{x,x), 



n(3) (x, y)^- -Giy, x)G{x, y)G{x, y) , 
6 

nW(x, y) ^'-^G{y, x)JV^J''zG^{x, z)G\z, y) . 

It is worth mentioning that the appearance of the gener- 
alized 5 function in the first local term is a consequence 
of the form of the effective action and the functional dif- 
ferentiation rule ([5]). For each vertex in the loop dia- 
grams there is a corresponding integral in the effective 



action. Because of the appearance of the generalized 5 
functions two of the integrals can be carried out trivially 
after functional differentiation. Further integrals persist 
in the self-energy. That is, four- and higher-loop contri- 
butions to n(x, y) contain integrations over space-time 
with the corresponding number of ^/—g factors to ensure 
the invariance of the self-energy. 



IV. KADANOFF-BAYM EQUATIONS 

Convolving the Schwinger-Dyson equations ([M)) with 
G from the right and using p^ we obtain 

i[U, + M^]G{x,y)^53i^x,y) 

+ Jy/^d-^zIl{x,z)G{z,y). (27) 

Next, we define the spectral function 

Gp{x,y)=i{Mx),^{y)]^), (28) 

and the statistical propagator 

GF{x,y)^^Mx),^{y)]+). (29) 

As is clear from the definitions, the statistical propagator 
of real scalar field is symmetric whereas the spectral func- 
tion is antisymmetric with respect to permutation of its 
arguments. For a real scalar field Gpix^y) and Gp{x,y) 
are real- valued functions 7]. The full Feynman propa- 
gator can be decomposed into a statistical and a spectral 
part 

G(x, y) = Gf{x, y) - §sign(x" - y")Gp(x, y) . (30) 

Upon use of the sign- and 5-function differentiation 
rules, the action of the D^; operator on the second term on 
the right-hand side of ([50)1 gives a product oig^^5{x^ ,y^) 
and S/qGp{x, y). Using the definition ([28)) and the canon- 
ical commutation relations in curved space-time [34I 



lim [(/9(x",x),7r(a;°,j7)]_ = iS{x,y), 



(31) 



where"* tt = g'^'^\/—gVoip, we find for the derivative of 
the spectral function 



VSGp(x,y) 



H^; y] 



(32) 



Multiplication of ([5^ by g^^5{x^ , y^) then gives the gen- 
eralized 5 function 6^{x^ y), which cancels the generalized 
6 function on the right-hand side of ([77)1 . 



* To simplify the calculation we set goi = 0. The off-diagonal 
components of the metric tensor can always be set to zero by an 
appropriate choice of the coordinate system ISal . Examples are 
the longitudinal and synchronous gauges. In the FRW universe 
this condition is fulfilled automatically. 



The local term of the self-energy ((26|) . proportional to IIf{x, y), in complete analogy to ([30l 
the S function, can be absorbed in the effective mass 



A 



Ar(x) =Ar + -G{x,x) 



(33) 



The remaining part of the self-energy can also be split 
into a spectral part, Ilp{x,y), and a statistical part, 

I 



Integrating along the closed time path in the direction 
indicated in Fig.[2l and taking into account that any point 
of the negative branch is considered as a later instant 
than any point of the positive branch, we finally obtain 
the system of Kadanoff-Baym equations: 



V X 

[D, + M\x)]GF{x,y) = fv^d^zUF{x,z)Gp{z,y)~ f ^^d^zllp{x, z)GF{z,y) , (34a) 



[D, + M%x)]Gp{x,y) = - f V^d^zUpix,z)Gpiz,y) . (34b) 



Comparing with the Kadanoff-Baym equations presented 
in [3, Q , we conclude that (p4|) appear to be the covari- 
ant generalization of the Kadanoff-Baym equations in 
Minkowski space-time. 

Equations p4p are exact equations for the quantum 
dynamical evolution of the statistical propagator and 
spectral function. It is important that, due to the char- 
acteristic memory integrals on the right-hand sides, the 
dynamics of the system depends on the history of its evo- 
lution [34J. 

To complete this section we derive explicit expres- 
sions for the spectral and statistical self-energies. Using 
symmetry (antisymmetry) of the spectral and statisti- 
cal propagators with respect to permutation of the argu- 
ments, we obtain for the three-loop contribution to the 
self-energy components: 



n?^(x,y) 



Tlf\x,y)^ 
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[Gpix, y)GF{x, y)GF{x, y) 



jGF{x,y)Gp{x,y)Gp{x,y)], (35a) 

--[3Gf{x, y)GF{x, y)Gp{x, y) 
\Gp{x,y)Gp{x,y)Gp{x,y)]. (35b) 



Four- and higher-loop contributions to the self-energy 
components contain integrations over space-time with x^ 
and y" as the integration limits. Introducing 



G4F(a;,2/) = /g y^d'^zGF{x,z)Gp{x,z) 

x[GUz,y)~lGl(z,y)] + {x^y}. 

G4p(x,j/) = /g y^d zGf{x,z)Gp{x,z) 
X [2Gf(2, y)Gp{z, y)] -{x^y} , 



(36a) 
(36b) 



and spectral components of the self-energy as 
n^ (a;,y) = — [GF{x,y)GiF{x,y) 



\Gp{x,y)Gip{x,y)], (37a) 



A3 



we can write the four-loop contribution to the statistical 



nf\x,y) ^ — [GF{x,y)Gip(x,y) 

+ Gp(x,y)GiF{x,y)]. (37b) 

Of course, all quantities entering the Kadanoff-Baym 
equations must be renormalized. The renormalization 
at finite temperature has been developed in [35|, |36|, |37|, 
|38| . A generalization to out-of-equilibrium systems with 
non-Gaussian initial conditions has been obtained in |39l . 
|40[. A renormalization procedure at tadpole order in 
the Gaussian scheme in the expanding universe has been 
applied to the analysis of Kadanoff-Baym equations in 



V. QUANTUM KINETICS 

Introducing the retarded and advanced propagators 

Gi,(a:,2/) = 0(x°-y°)Gp(x,2/), (38a) 

G^(x,2/) = -0(yO-x")Gp(x,y), (38b) 

and the corresponding definitions for the self-energies, 
one can rewrite the system of Kadanoff-Baym equations 
in the form: 

[U, + M^{x)]GF{x,y) = -JV^d*z9{z°) 

X [Uf{x, z)Ga{z, y) + TIbIx, z)Gf{z, y)] , (39a) 

[a^ + M\x)]Gp{x,y) = -Jy^d^zeiz") 

X [np{x, z)Ga{z, y) + nfl(x, z)Gp{z, y)] . (39b) 



The system (|39|) should be supplemented by the analo- 
gous equations for the retarded (advanced) propagators; 
they can be derived from (|34b|l upon use of ([32|) 

[n, + M\x)\GR^A){x.v)^53{x,y) 

-JV^d^znR^A){x,z)GR(^A)iz,y)- (40) 



Let us now interchange x and y on both sides of the 
Kadanoff-Baym equations ([5^ . The difference (sum) of 
the original and resulting equations are referred to as 
the kinetic (constraint) equations for the spectral func- 
tion and the statistical propagator. Using the relation 
G_R,(x,2/) — GA{y,x) and symmetry (antisymmetry) of 
the statistical propagator (spectral function) we obtain 



[D, T Oy + M\x) T M\y)]GF{x, y) = 

= -JV^d'^ze{z")[nF{x,z)GA{z,y) T GR{x,z)nF{z,y) 

+ Unix, z)Gf{z, y) T Gf{x, z)T\a{z, y)] , 
[D, T Dy + M\x) T M\y)]Gp{x,y) = 

= -/x/^d*z^(z")[np(x, z)Ga{z, y) T Gr{x, z)Iip{z, y) 

+ ^b{x, z)Gp{z, y) =F Gp(x, 2)0^(2, y)] , 



(41a) 



(41b) 



Interchanging x and y on both sides of the equation for Ga{x, y) and adding it to the equation for Gf!.{x, y) we obtain 
the constraint equation for the retarded propagator: 



p, + ny + M^{x) + M^{y)]GR{x, y) 



253 {x, y) - S./^d^z[nR{x, z)Gr{z, y) + Gr{x, z)\Ir{z, y)] 
I 



(42) 



Next, we introduce center and relative coordinates. In 
Minkowski space-time they are given by half of the sum 
and by the difference of x and y, respectively [9|. In 
other words the center coordinate lies in the middle of 
the geodesic connecting x and y, whereas the relative 
coordinate gives the length of the "curve" ^ connecting 
the two points. 

Consider now curved space-time. Let <j be the affine 




FIG. 3: Arrangement of the points along the geodesic. 

parameter of the geodesic connecting x and y (see Fig. [3]) 
and ^(<;) a function mapping ? onto the points of the 
geodesic, with 



X 



iW), y" = e(0. 



(43) 



The center coordinate lies in the middle of the geodesic, 
i.e. it corresponds to <^x = ^W + "?")• The relative 
coordinate is given by the sum of the infinitesimal dis- 
tance vectors df " along the geodesic, all of which must 



have been submitted to parallel transfer to ?x from the 
integration point on the curve^. According to [431 this 
implies 

(44) 

All quantities in equations (|4ip are now recast in terms 
of X" and s". Up to higher order, proportional to the 
curvature tensor terms, the Laplace-Beltrami operator is 
given by [4a | 



□. 



-D"D^ 



'' ±D'>' 



4 " ds'^dsa dsa ' 

where Da is the covariant derivative 



(45) 



(46) 



In Minkowski space-time geodesies are straight lines. 



Calzetta and Hu Il7ll4l| have employed a different method based 
on the use of Riemann normal coordinates and the momentum 
representation of the propagators. Their approach has some ad- 
vantages for the study of the quantum kinetics equations. Here 
we are mainly interested in the Kadanoff-Baym and Boltzmann 
equations and consider the derivation of the quantum kinetic 
equations as an intermediate step connecting both of them. For 
this reason, we adopt the covariant definitions of the midpoint 
and distance vectors introduced by Winter [4y|, which allow us 
to keep the analysis manifestly covariant in every step. 



Note that in (|45p we have neglected the corrections pro- 
portional to the Riemann and Ricci tensors. Next, we 
Taylor expand the effective masses to first order around 
the center coordinate X 



M^ 



M^{X)±y^DaM^{X) 



(47) 



where the minus sign corresponds to y whereas the plus 
sign corresponds to x. The propagators on the left-hand 
side of (J4T|) can also be reparameterized in terms of the 
center and relative coordinates: GF{x,y) -^ Gf{X,s) 
and Gp{x,y) -^ Gp{X,s). 

On the right-hand sides we have convolutions of func- 
tions of X and z and functions of z and y. That is, we 
have to introduce the corresponding center and relative 
coordinates and perform the integration. Making use of 
the identity 

w + n = (<?' + <?") + (^* - <?") = 2ix + i^*~ <?") 

and Taylor expanding around <jx , we obtain to first order 

IIf{x, z) = IIf{X^^, S^cz) ~ n_F(X, Sxz) 



8^°' di; du" dq 



(48) 



Using furthermore the definition of the four-velocity and 
the geodesic equation 



we can rewrite 



in the form 



-viyu'' , 



(49) 



^f{x,z) K.IVf{X,S^:, 



2°zy 



dJVf{X,s,,), (50) 



where s"^ = (<;* — q")u°'{<,x)- Making use of the identity 

(^" + ',*) = W + ?") - (^' - ^*) = 2 ^x - W - C) 

we get a similar expression for the functions of z and y 

Ga{z, y) « Ga{X, s,y) - \sZ,D^Ga{X, s,y) . (51) 

To perform the integration of the product of ([50]) and 
((5T|) , we shift the coordinate origin to <;x and replace the 
integration with respect to z by integration with respect 
to distance sxz from X to z along the geodesic. More- 
over, we approximate^ \F~9z ^^ ^^^ value at the origin 



'' The next-to-leading term of the Taylor expansion is propor- 
tional to the convolution of the Christoffel symbol 3l^j , \/— 3^ ~ 
■^— g^(l + r^^s"). This correction can in principle be taken 
into account and would induce additional terms proportional to 
id /dp" on the right-hand side of the quantum kinetic equation. 
Since such term are neglected in the Boltzmann approximation, 
the collision terms do not receive any corrections. 



The Kadanoff-Baym equations describe the dynamics 
of a system in terms of the spectral function and statis- 
tical propagator. The latter ones are functions of two 
coordinates in the four-dimensional space-time. By in- 
troducing center and relative coordinates we have traded 
one set of coordinates for another one. Performing the 
so-called Wigner transformation, one can also trade one 
of the arguments defined in the coordinate space for an 
argument defined in the momentum space. In curved 
space-time [42| 

GF{X,p) = ^^xJd^se'P'GF{X,s), (52a) 



GFiX,s) = 



9xJ (2^) 



e-'P'GF{X,p). (52b) 



Note that in ([52)) and in the rest of the paper we use 
contravariant components of the space-time coordinates 
and covariant components of the momenta. Let us also 
note that 



dm 



1 



p 



~9x (2^)' 



is the invariant volume element in momentum space. The 
definition of the Wigner transform of Gp{X,s) differs 
from (|52ap by a factor of —i so that Gp{X,p) is again 
real valued. 

As follows from (j52bp . differentiation with respect to 
s" is replaced after the Wigner transformation by p"' 



_d_ 
9s° 



-ip^ 



(53) 



Upon integration by parts we also see that s" is replaced 
by differentiation with respect to p" : 



. d 
'dp^ 



(54) 



Consequently the Wigner transformed covariant deriva- 
tive reads 

i^.-^.-^ + rf,p,A. (55) 

Correlations between earlier and later times are expo- 
nentially suppressed, which leads to a gradual loss of the 
dependence on the initial conditions [9|, |3J| . Exploiting 
this fact, one can drop the 9 function from the integrals 
in the difference equations (|¥l]) . Furthermore we let the 
relative-time coordinate s*^ range from — cx) to oo in or- 
der to perform the Wigner transformation, see [3J, |43| 
for a detailed discussion of these approximations. Then 
using ([51)1 and ([55)1 we obtain for the Wigner transform 
of the first term on the right-hand side of (|41ap : 

IV-9zd'^z HfIx, z)Ga{z, y) -^ 
nF{X,p)GA{X,p) + i^{nF{X,p),GA{X,p)}pB , (56) 
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where the Poisson brackets are defined by 
{AiX,p),B{X,p)}pB = ■^A{X,p)V^BiX,p) 

~Va.AiX,p)^BiX,p). (57) 

OPa 

Comparing (j57p to its Minkowski-space counterpart we 
see that the derivatives with respect to X are replaced 
by the covariant derivatives, just as one would expect. 

Wigner transforming the rest of the terms we obtain a 
rather lengthy expression which can be substantially sim- 
plified with the help of the relations between Gii{X,p), 
Ga{X,p), and Gp{X,p). Recalling the Fourier transform 
of the 6 function, 

fds°exp(iws°)0(±s°) = lim ^^, 



we find that 



GRiX,p) = - 
Ga{X,p) = - 



dhjGp{X,p,uj) 
2t: po — Lu + ie 
dhJ Gp{X,p, Lo) 



(58a) 
(58b) 



2tt pq — UJ — ie 
From comparison of (j58bp and (|58ap it follows that 

GAiX,p) = GUX,p). (59) 

Recalling furthermore that the S function can be approx- 
imated by 



6(u;) — lim , „ ^ , 

^ ^ e^+o 7r(a;2 + ^2) 



(60) 



we also find that 

GR{X,p)-GAiX,p) = tGp{X,p). (61) 

Analogous relations also hold for the retarded and ad- 
vanced components of the self-energy. 

As can be inferred from (|45)) and (|47p . the Wigner 
transform of the left-hand side of (HTI) reads ^ 



D,-ny + M^{x)-M^iy) 



i 2p'^Va + DaAf' 



d 

dpa 



(62) 



Introducing the quantity 

niX,p)=p^Pp~M^{X)-I[h{X,p), (63) 



Additional contributions arising from the decomposition of the 
Laplace-Beltrami operator are proportional to Riemann and 
Ricci tensors and to the curvature (see Eq. (4.40) in 1431 ) and 
may be relevant in strong gravitational fields. Since all these 
terms contain at least one id/dpa derivative, they do not con- 
tribute in the Boltzmann approximation. 



where Iih{X,p) = Renfl(A', p), and collecting the terms 
on the right-hand side of the kinetic equation (PT|) one 
can write the kinetic equation for the Wigner transform 
of the statistical propagator in the compact form: 

{h[X,p),GF{X,p)}pB 

= GF{x,p)np{x,p) - nF{x,p)Gp[x,p) 

+ {nf(X,p),G,.(X,p)}pB, (64) 

where Gh{X,p) = KeGii{X,p). The same procedure 
leads also to a kinetic equation for the Wigner transform 
of the spectral function 



{nix,p),Gpix,p)} 



PB 



= {np{x,p),Gh{x,p)}pB. 



(65) 



As has been mentioned in the previous section, the ex- 
act quantum dynamical evolution of the system depends 
on its whole evolution history. Mathematically, this man- 
ifests itself in the memory integrals on the right-hand 
sides of ((34|) . In fact, performing the linear order Taylor 
expansion around X, we take into account only a very 
short part of the history of the evolution. Since the ex- 
pansion coefficients are defined at X, after the integration 
we obtain equations which are local in time. 

Next we consider the Wigner transform of the con- 
straint equation for the retarded propagator (|^ . On the 
left-hand side we have D^; -I- Dj, = 2dsa dg^ , to first order 
in the covariant derivative, whereas M^{x) + M'^(y) « 
A'P(X). On the right-hand side the Poisson brack- 
ets cancel out and only the product of Ylji{X,p) and 
Gii{X,p) remains. Finally, the Wigner transform of the 
generalized S function is just unity. Therefore, we get 
an algebraic equation for the Wigner transform of the 
retarded propagator 



[P>M 



M^{X)-I[R{X,p)]GRiX,p) 



4. 



(66) 



Equation ([66]) implies that the real part of the retarded 
propagator is given by 



Gh{X,p) 



-n{x,p) 



n^{x,p) + Wp{x,p) 



(67) 



Note that Gh{X,p) vanishes on the mass shell, which is 
defined by the condition Q.{X,p) = 0. As follows from 
([55)1 and (pT|) , the Wigner transform of the spectral func- 
tion is twice the imaginary part of the retarded propaga- 
tor: 



Gp{X,p) 



-np(x,p) 



nHx,p) + \Til{x,p) 



(68) 



Equation (|68p is also a solution of (|65p . To first order 
in the covariant derivative the Wigner-transform of the 
constraint equation for the statistical propagator reads 

n{X,p)GF{X,p) = \{tiF{X,p),Gp{X,p)}pB 

+ \{GF{X,p),\lp{X,p)}pB+\lF{X,p)Gh{X,p). (69) 



The constraint equation for Gp{X,p) is no longer alge- 
braic and can not be solved analytically in general. How- 
ever, let us assume for a moment that the system is in 
thermal equilibrium. In this case all the quantities are 
constant in time and space and the Poisson brackets in 
([69|) vanish identically. The solution of the resulting al- 
gebraic equation then reads 



G7ip) 






G7ip). 



(70) 



That is, we have obtained the fluctuation-dissipation re- 
lation. It only remains to calculate the ratio of the spec- 
tral and statistical self-energies. This can be done us- 
ing the relation (PD|) and the KMS periodicity condition, 
G{x,y)\x=o = G{x,y)\x=-i(3, where /? is the inverse tem- 
perature. Wigner-transforming this equation and using 
([70)1 we obtain 



Gp\p) 



Aeq) 



(P)- 



&:'^\p). 



(71) 



where rv-'^'^' is the Bose-Einstein distribution function. 

To complete this section, we have to express the 
Wigner transforms of the spectral and statistical self- 
energies in terms of the Wigner transforms of the spectral 
function and statistical propagator. Using the definitions 
of the Wigner transformation and its inverse we find for 
the Wigner transform of a product of functions of the 
same arguments: 

fi{x,y)---fn{x,y) -» fi...fn{X,p) 

= JdU^^ . . . dnlJ27rfV^^S^{-p + p, + ...p„) 

x/(X,pi).../(X,p„). (72) 

Note that 5g{q) = \f^^x ^i'i) represents the momentum- 
space generalization of the 6 function, invariant under co- 
ordinate transformations (this can be checked with help 
of the scaling property of the S function). Keeping in 
mind that the definition of Gp{X,p) contains an addi- 
tional factor of —i we can then write the Wigner trans- 
forms of (1551) in the form 



n'^^\x,p) 



A 



^,GUX,p) + IGfGI{X,p)], (73a) 



ufHx,p) = - -[3G%Gp{X,p) + iG3(X,p)] . (73b) 

The expression for the Wigner transform of the three- 
loop retarded self-energy can be obtained from ([75]) by 
replacing one of the Gp by Gr. The Wigner transforms 
of the four-loop contributions (j35p can be written in a 
similar way 

n^^\x,p) ^^[G^p{X,p) + \G^p{X,p)] , (74a) 
^i^Hx,p) ^^[G^4p{X,p) + G^f{X,p)] . (74b) 



Note, however, that Gaf and G^p are Wigner transforms 
of convolutions of four two-point functions, 

GiF{x,y) = J^^d*zGFix,z)GR{x,z) 

x[Gl{z,y)^lGliz,y)] + {x^y}, (75a) 
Gip{x,y) ^ J^/^d'^zGFix,z)GR{x,z) 

X [2Gf{z, y)Gp{z, y)] - {x ^ y} , (75b) 

where we have used the definitions of the retarded and 
advanced propagators and dropped again the 9{z°) fac- 
tor. Proceeding as in Eq. (|56p and making use of the 
relations ([59|) and ([6T|l . we obtain for the Wigner trans- 
forms of GiF and G^p 

G,f{X,p) = 2[G|(X,p) -I- \Gl{X,p)]G7m^R]{X,p) 
+ i{G|(X,p) + \Gl{X,p),ChMp{X,p)}p,B. , (76a) 



Gip{X,p) = 4GfGp(X,p)Gf5R[Gk](X,p) . 



(76b) 



Finally, the expression for the Wigner transform of the 
four-loop retarded self-energy can be obtained from 
(|74bp by replacing Gp with Gr and G/^p with G4^r. The 
latter one is related to G^p by Eq. (|58ap . 



VI. BOLTZMANN KINETICS 



The spectral function ([55]) has approximately Breit- 
Wigner shape with a width proportional to the spectral 
self-energy. The area under Gp{X,p) is determined by 
the normalization condition, 



^00 



27r 



Gp{X,p)pQdpo = 1, 



(77) 



which is a direct consequence of (1321) and the antisymme- 
try of the spectral function with respect to permutation 
of its arguments. In the limit of vanishing coupling con- 
stant the width of the spectral function approaches zero, 
whereas its on-shell value goes to infinity, see Eq. (|68p . 
Equation (j60p then implies that in this limit the spectral 
function takes the quasiparticle form p| 



Gp{X,p) = 27rsign(po) S {g^^y^p, - W 



(78) 



Note that ([75]) is consistent with the normalization condi- 
tion ((77)) . The signum-function appears in ((78)) because 
np(X, p) is an odd function of po- Since the magnitudes 
of lip, n^ and of the local term of the self-energy are 
controlled by the same coupling we have also neglected 
them in ^(X^p). In the same limit Eq. ((S5)) for the spec- 
tral function simplifies to 



p"I?„Gp(X,p) = 



(79) 



and indeed admits a quasiparticle solution ((75| . Note 
that Eqs. ((79)) and ((75)) state that the effective mass M 
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of the field quanta does not change as they move along 
the geodesic, just like it is the case for particles. 

Motivated by the fluctuation-dissipation relation (fflj) 
we can trade the statistical propagator for some other 
function: 



GFiX,p)^ [niX,p) + ^]GpiX,p) 



(80) 



However, if both Gf{X,p) and Gp{X,p) are smooth func- 
tions then relation (|80p is merely a definition of n(X,p). 
In the quasiparticle approximation the spectral function 
is divergent and forces the momentum argument of n to 
be on the mass shell. For this reason the quasiparticle ap- 
proximation for the statistical propagator (fSO ) is usually 
referred to as the Kadanoff-Baym Ansatz [a, l3ll | . 

Let us now tentatively put the coupling constant to 
zero. In this case the right-hand sides of the kinetic 
equations ((64|) and ((65|l vanish. In this case Gf{X,p) 
and Gp {X, p) are constant in space and time even if the 
system is out of equilibrium. If we now "increase" the 
coupling constant again, then the gain and loss terms 
on the right-hand side of (|64p will induce nontrivial dy- 
namics for the statistical propagator. This in turn will 
induce a time and space dependence of the spectral and 
statistical propagators thus leading to nonvanishing Pois- 
son brackets on the right-hand sides of the kinetic equa- 
tions. The magnitude of the derivatives with respect to 
the time and space coordinates are therefore proportional 
to some (positive) power of the coupling constant. Con- 
sequently the contribution of the Poisson brackets in ([M]) 
is effectively of higher order in A than the contribution 
of the gain and loss terms. These considerations justify 
the dropping of the Poisson brackets and of the local and 
nonlocal contributions to the effective field mass in the 
kinetic equations. In other words, they legitimate the use 
of the quasiparticle approximation.^ 

As has been argued above, the Poisson brackets par- 
tially take into account the memory effects. Neglecting 
the Poisson brackets we completely ignore the previous 
evolution of the system. Physically this corresponds to 
the Stosszahlansatz of Boltzmann. 

From Eqs. dMl), GH) and (EHl) it follows that in this ap- 
proximation the kinetic equation for the statistical prop- 
agator turns into an equation for the evolution of the 
one-particle distribution function n{X,p): 

[p"I?„n(X,p)]Gp(X,p) 

= i [n> (X, p)G< (X, p) - G> (X, p)n< (X, p)] , (81) 

where we have introduced 

G^{X,p) = Gf{X,p)±\Gp{X,p) (82) 



® If we were interested in higher order processes, for instance in 
the 2 — > 4 scattering which is of the fourth order in the coupling 
constant, we would have to use the so called extended quasipar- 
ticle approximation [ij, HI, HI, 113, El, H^l . Since in this paper 
we limit ourselves to the processes of at most third order in A, 
the quasiparticle approximation is sufficient for our purposes. 



and their self-energy analogs n^(X,p). The symmetry 
(antisymmetry) of the statistical (spectral) propagator 
with respect to permutation of its arguments and the 
definition of the Wigner transformation imply that 

Gf{X,p) = GpiX, -p), Gp{X,p) = -Gp{X, ~p) . (83) 

Therefore, for a single real scalar field, we have 

G^(X,-p)-%(X,p), (84) 

and a similar relation for the self-energies. 

Explicit expressions for n> (X, p) can be obtained after 
some algebra from Eqs. ([75)1 and ([7^ . For illustration 
purposes we first derive n> (x, y) and then perform the 
Wigner transformation. Using the decomposition 

G(x, y) - 0(xO - yO)G> (x, y) + 0(y" - x^)G< {x, y) (85) 

we obtain for the three-loop contribution 



t(3)^ 



A2 



T]}^'{x,y) = -—G^{x,y)G^{x,y)G^{x,y). (86) 
Its Wigner transform reads 

nf{x,p) = - ^Sdnldn\<mt(2i:)Hg{-p -t + k + q) 



X G<(X,t)G>(X,fc)G>(X,q) 



(87) 



where we have used relation (|84|) . It describes 2 <-> 2 
scattering and corresponds to the tree-level Feynman di- 
agram in Fig.m 

k P 




q t 




FIG. 4: Feynman diagrams of 2 
one-loop levels. 



2 scattering at tree and 



Expression for the four-loop contribution contains in- 
tegration over the contour 

A3 



tif{^,y) = yG^(x,2/)/V^,d4z0(zO)x 

[Gf{x, z)Gr{x, z)G|(2;, y) + G| (x, z)Ga{z, y)GF{z, y)]. 

After some algebra we obtain for the Wigner transform 
of (|88p in the Boltzmann approximation (that is, with 
the Poisson brackets neglected) 

A3 



uf{x,p) = '-jdutdn^^dnU2n)\i~p- t + k + q) 

X % {X, t)G^ (X, k)G^ (X, q) L{X,k + q), (89) 



where 

L{X,p) = Jdntdn^^i2n)^Sg{-p + k + q) 

y.2GF{X,k)Gh{X,q). (90) 

From ([55)1 it follows that L{X,p) is the same for the 
forward and inverse processes. As is demonstrated in 
Appendix El it corresponds to the integrals of the one- 
loop Feynman diagrams in Fig.[51 

Let us note here that the contribution(s) of a partic- 
ular term of the 2PI effective action to the Boltzmann 
equation can be deduced by cutting the 2PI diagrams by 
a connected line in all possible ways. The three-loop con- 
tribution, for instance, can be cut in only one way and 
the result can be represented as a product of two tree- 
level scattering diagrams. The four-loop contribution 
can be cut in three equivalent ways and the result can 
be represented as a product of tree-level and one-loop 
scattering diagrams, see Fig. [5] There are two five-loop 






FIG. 5: The correspondence between the diagrams contribut- 
ing to the 2PI effective action and the contributions to the 
Boltzmann collision terms. 

loop contributions to the effective action [3l|. Apply- 
ing the same procedure to one of them we would obtain 
interference terms of two one-loop scattering diagrams 
and interference of tree-level and two-loop scattering di- 
agrams. Cutting the second, "eye" , diagram we would 
obtain interference of tree-level and two-loop scattering 
diagrams and also interference of two 2 — > 4 diagrams. 

The quasiparticle approximation (|78p for GpiXjk) in 
([50)1 forces one of the intermediate states in the loop to be 
on the mass shell. On the contrary Gh, which describes 
the second intermediate state in the loop, vanishes on 
the mass shell. That is, the real intermediate state con- 
tributions (2^2 scattering into two on-shell states fol- 



J 
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lowed by another 2 — > 2 scattering) are automatically 
subtracted from the four-loop self-energies. 

Also note, that initial and final states and on-shell 
intermediate states can be clearly distinguished in this 
formalism: the former ones are described by G^ compo- 
nents, whereas the latter ones by Gp or Gp components. 
Performing the integration and taking into account that 
one of the intermediate states is on-shell, we obtain the 
following expression for the loop integral: 



L{X,p) = lim/ 



dk 2n(A:,k) + l 



0-^ (27r)3 
p^ — 2pk 



2Ek 



p^ + 2pk 



(p2_2pfc)2 + e2 (p2 _^ 2pfc)2 + £2 



(91) 



where k = [Ek , k) is the on-shell four-momentum ex- 
pressed in terms of the "physical" components: Ek = 
ko/^/goo, etc. In ([?T|l the background plasma "affects" 
only one of the internal lines; the other one is off-shell 
and we can not associate the particle number density 
with it. 

Next, we integrate the left- and right-hand side of ([8T1) 
over po and choose the positive energy solution of fTS]) 
on the left-hand side. On the right-hand side both, the 
positive and the negative energy, solutions contribute. 
For positive Pq momentum-energy conservation allows 
the following three combinations: 

a) fco > 0, 9o > 0, to > , 

b) fco > 0, go < 0, to < , 

c) fco < 0, 9o > 0, to < . 

As far as the three-loop self-energy ([57)1 is concerned, 
each combination leads to the same result, i.e. an over- 
all factor of 3 appears. For the four-loop self-energy 
the arising terms are not equal due to the presence of 
the loop integral L in ([50)1 . Taking this into account and 
comparing ([87|l and ([89]) we see that in the 2PI formalism 
the effective coupling at nonzero particle number density 
at one-loop level contains a sum of three L{X,p) func- 
tions with the arguments corresponding to s-, t— and 
M-channel scattering: 



A'^{X,k,q,t) 



; A2(l-A[L(A,fc + q) 
L{X,k-t)+L{X,q- 



t)] 



(92) 



After some algebra, the use of (|84p and redefinition of 
the momenta we finally arrive at the Boltzmann equation 
for the distribution function: 



p"I?an(A,p)=-^y 



rfk dq 



dt 



(2^)3 (2^)3 EkE.Et 



S{Ep + Et-Eq~ Ek)5{p + t - q - k) v42(X,k,q,t) 



X {n{X, p)n(X, t) [n{X, k) + 1] [n{X, q) + 1] - [n{X, p) + 1] [n{X, t) + l]n{X, k)n(X, q)} 



(93) 
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It is interesting, that the only remnant of the curved 
structure of space-time is the covariant derivative on the 
left-hand side of the Boltzmann equation. In the case of 
greatest practical interest - the Friedmann-Robertson- 
Walker universe - it takes the form 



p"X>„n 



E f d 



a \drj 



i^U'- «4' <-) 



where 77 is the conformal time. An integral form of the 
Boltzmann equation in the FRW universe as well as in a 
space-time with linearly perturbed FRW metric can be 
found, for instance, in [50|. 

On the right-hand side, all the y/—gx factors have 
disappeared due to the introduction of the "physical" 
momenta and energies. In other words, the transition 
amplitudes in the scattering terms are independent of 
the space-time metric, which justifies many earlier cal- 
culations. It is also remarkable that if only pointlike in- 
teractions (i.e. only the three-loop contribution to the 
2PI effective action in the considered case) are taken into 
account, Eq. (|93p coincides with the classical Boltzmann 
equation with the collision term calculated in vacuum. 
The inclusion of four- (and higher-loop) corrections to 
the effective potential induces further terms in the Boltz- 
mann equation. These terms correspond to the remnant 
space-time integrals in the self-energy and involve ad- 
ditional momentum integrals over the distribution func- 
tions. 



VII. SUMMARY AND CONCLUSIONS 

In this paper we have considered the dynamics of 
an out-of-equilibrium quantum system in a background 
gravitational field in the Schwinger-Keldysh formalism . 
As one would expect, the resulting equations turned out 
to be covariant generalizations of their Minkowski-space 
counterparts. 

Remarkably, in the Boltzmann approximation the only 
remnant of the curved structure of the space-time is the 
covariant derivative on its left-hand side. The matrix 
elements of the scattering terms on the right-hand side 
are independent of the metric. This justifies earlier cal- 
culations where this has been assumed implicitly. Fur- 
thermore, if only the tree-level processes are taken into 
account, then the resulting equation coincides with the 
Boltzmann equation with the collision term calculated in 
vacuum. Processes described by loop diagrams, which in- 
duce corrections to the self-coupling, involve additional 
momentum integrals over the distribution functions, so 
that the resulting contributions differ from those calcu- 
lated in vacuum. 

Interestingly, loop corrections, i.e. processes with in- 
termediate off-shell states, can be taken into account 
even if the quasiparticle Ansatz is applied. As far as on- 
shell intermediate states are concerned, there is a clear 
distinction between them and the initial and final states: 
the former ones are described by Gp (or Gp) components, 



whereas the latter ones are given by G^ components, see 
Eqs. ((90| and (l93l) . It is important that in the used for- 
malism the problem of double-counting, which is cured 
by a real intermediate state subtraction procedure in the 
standard approach, does not arise at all. 

For leptogenesis, this implies that whereas the washout 
processes described by contact interactions (they are 
present for instance in the supersymmetric extensions of 
the Standard Model) can be treated essentially classi- 
cally, the correct treatment of the decay processes (which 
generate the asymmetry) and the scattering processes 
mediated by the right-handed neutrino (which washout 
the asymmetry) requires the use of the Kadanoff-Baym 
approach. 

Since the peculiarities of the calculation, related to the 
presence of a background gravitational field, are deter- 
mined only by transformation properties of the fields - 
scalar fields in the present case - the developed formal- 
ism can be applied to arbitrar y sy stems of scalar fields 
without any modifications. In [5ll | we study further im- 
plications of this formalism for leptogenesis and calculate 
the vertex contribution to the CP-violating parameter 
at nonzero particle densities in the framework of a toy 
model that qualitatively reproduces the features of pop- 
ular leptogenesis models. The analysis of the self-energy 
contributions to the CP-violating parameter will be per- 
formed in I52I. 
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APPENDIX A: 2 ^ 2 SCATTERING 

The tree-level amplitude of 2 ^^ 2 scattering (see 
Fig.[¥]) in Minkowski space-time is given by 



j^jtree ^ _^^ 



(Al) 



There are also three one-loop diagrams which contribute 
to the scattering amplitude. Their contribution reads 



M 



loop 
fi 



-A2 



d'^Cd^riSi-a + ^ + T]) 



2(27r)4 J [^2 _ ]^.p _^ jg][^2 _ ]^.p j^ ^ 



(A2) 



where a is equal to k + q, to k — t or q — t (see Fig.H]). 
Because of the presence of the (5-function one of the inte- 
grations (for instance, over ry) can be performed trivially. 
Calculating residues of the integrand we can perform the 
integration over d^Q. The result of the integration reads 



m: 



loop 



iX' 



2(27r) 



d^^ 
2EI 



1 



1 



$2 + 2^(7 C'-2CfT 



(A3) 
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The quantity which enters the right-hand side of the 
Boltzmann equation is the ampHtude modulo squared. 
To leading order in small A it is given by 

+ L™'=(fc-t) + L™=(g-i)]), (A4) 

where L"''^(cr) coincides with (PTjl if n{X, k) and e are set 
to zero. The former condition arises from the fact that in 



this Appendix we calculate the scattering amplitudes in 
vacuum, whereas the latter one is related to the fact that 
we have not subtracted the contributions of real interme- 
diate states to the one-loop amplitude. Comparing (|A4p 
with (|92p we conclude that L{X,p) indeed describes the 
integrals of the one-loop diagrams. 
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